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O ' Abstract 
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Graph-theoretical approach is used to study cluster formation in mesocsopic systems. 
Appearance of these clusters are due to discrete resonances which are presented in the 
form of a multigraph with labeled edges. This presentation allows to construct all non- 
isomorphic clusters in a finite spectral domain and generate corresponding dynamical 
systems automatically. Results of MATHEMATICA implementation are given and two 
possible mechanisms of cluster destroying are discussed. 
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1 Introduction 

> ' 

Q\ \ Mesoscopic regimes are at the frontier between classical (single waves/particles) and statistical 
^ [ (infinite number of waves/particles) description of physical systems. Mesoscopic systems is 
^ ■ very popular topic in various areas of modern physics and can be met in wave turbulent theory, 
O ■ condensed matter (quantum dots), sociology (opinion formation), medicine (dynamics of cardio- 

O 



vascular system), etc. For instance, statistical wave turbulence theory is based on Kolmogorov's 
suggestion on spatial evenness of turbulence and does not describe observed organized struc- 
tures extending over many scales like boulders in a waterfall. Also many laboratory experiments 
stay unexplained in terms of statistical approach as in pQ where the experimental results have 
^ ■ been presented for water turbulence excited by piston-like programmed wave-makers in water 
^ . flume with dimensions 6 x 12 x 1.5 meters. The main goal of this experiments was to establish 
^ ! a power-law scaling for the energy spectrum, E ~ 0~ u , with some fixed v coming from statis- 
tical considerations and being wave dispersion function. It turned out that discrete effects 
are major and statistical predictions are never achieved: with increasing wave intensity the 
nonlinearity becomes strong before the system loses sensitivity to the discreteness of spectral 
space. 

Our next example is taken from a quite different area of research - sociology. Finite size 
effects in the dynamics of opinion formation have been profoundly studied quite recently in [2] 
with essentially the same conclusions made. Namely, some changes of a system can be observed 
only when "finite number of agents in the model takes a finite value" and thermodynamic limit 
does not describe behavior of these systems. It was shown that resonance by which a finite-size 
system is optimally amplified by a weak forcing signal (identified as an advertising agent) is 
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determined by the size of mesoscopic system and the largest peak of the spectral density is 
observed at the driving frequency. 

Another interesting example of a mesoscopic system can be found in [3] where the flow 
of blood through the system of closed tubes — the blood vessels - - is described by wave 
equations. A model of the cardiovascular system as a system of coupled oscillators is proposed 
and conditions of their resonance are studied. 

Now the question is: what have boulders in waterfalls, advertisement and cardiovascular 
system in common? From mathematical point of view the answer is very simple: dynamics of 
all these systems can be interpreted as discrete resonances. Notice that resonance conditions 
have the same general form for wave and quantum systems (see, for instance, [I] for 4-photon 
processes); and have to be studied in integers. In this paper we have chosen a wave turbulent 
system as our main example and therefore use wave terminology. 

From now on we regard resonance conditions of the form 

01±0 2 ±...±0 S = O, k 1 ±k 2 ±....±k s = Q (1) 

where 0i = 0(h), s < oo, with k and 0i = 0{ki) being correspondingly wave vector and 
dispersion function. Specific features of these systems described by Fourier harmonics with 
integer mode numbers were first studied in [5] (we call them further discrete wave systems, 
DWS). It is well-known that fully statistical description of a wave turbulent system yields wave 
kinetic equation [6] analogous to kinetic equation known in quantum mechanics. A counter 
part to kinetic equation in DWS is a set of few independent dynamical systems of ODEs on the 
amplitudes of interacting waves. The theory presented in [5] was based on a collection of pure 
existence theorems [7] and has been developed with the understanding that discrete effects are 
only important for small \k\ of order ~ 10 while in larger spectral domains statistical regimes 
do occur. Numerous results of the last few years ([8], [9], \±0\, pj], [12], [T3j, [TJ] just to 
mention a few of them) showed that the general conception - discrete effects are only important 
in small spectral domains - should be revised because these effects are in fact observed in the 
systems where thousands of Fourier harmonics are taken into account, i.e. in a wide range 
of mesoscopic wave systems. A model of laminated wave turbulence has been presented in 
[15] which explains the appearance of coherent structures in arbitrary big but finite spectral 
domains. This model put a novel computational problem of solving (CQ) in integers of order 10 3 
and more. Fast generic algorithms for a big class of irrational and rational dispersion functions 
have been developed in [TB] which can be used for a wide range of dispersion functions; they 
were also implemented (for cases s = 3 and s = 4) and computation time on a Pentium-3 is 
of order 5 to 15 minutes in computation domain \m\, \n\ < 10 3 , where k = (m,n), m, n 6 Z. 
Straightforward computations for the same examples (without using our algorithms) takes few 
days with similar computer and computation domains of order 10 2 . 

In this paper we study the structure of the solution set of ([T]) using graph-theoretical ap- 
proach and develop a special technique to construct all independent clusters and corresponding 
dynamical systems. We present the whole solution set as a multigraph with labeled edges so 
that each connected component of this multigraph correspond to a special dynamical system 
of ODEs on the wave amplitudes. The most important fact about is this construction is follow- 
ing: is provides simultaneous isomorphism of multigraph components and dynamical systems. 
Using algorithms [16] we have developed a MATHEMATICA program package (at present only 
for s = 3 and 2D-waves) capable to 1) construct all independent clusters of the solution set 
of ([I]) in a given computation domain; 2) draw them as a multigraph on a plane; 3) write out 
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explicitly all dynamical systems appearing in the chosen spectral domain. Some results of our 
implementation are given, possible directions for further research are briefly discussed. 



2 Discrete 3-wave resonances 

As our main example, 3-wave resonances covered by barotropic vorticity equation (BVE) has 
been chosen. This equation, also known as Obukhov-Charney or Hasegawa-Mima equation, is 
important in many physical applications - from geophysics to astrophysics to plasma physics: 
the equation was again and again re-discovered by specialists in very different branches of 
physics. In particularly, this equation describes ocean planetary waves 

with non-flow boundary conditions in a rectangular domain 

if) = for x G [0,LJ, y G [0,L y ], 

where (5 is a constant called Rossby number and < e <C 1 is a small parameter. A linear 
wave has then form [17] 

/3 urn nn 

A cos(— x + 0t) sin — — x sin — y, m,n G Z 
20 L x Ly 

and dispersion function can be written as 



= 2//9^)> + Q». 

After obvious re-normalization we write out resonance conditions for 3-wave interactions as 
follows: 

' - 1 I 1 - 1 

\fml+n\ yjm\+n\ \/ m 3+ n 3 (3) 

n\ ±n 2 = n 3 

This system will be our main subject to study in this paper. 

Just for completeness of presentation we present here a simple idea underlying our algorithm 
for computing integer solutions of (for more details see [16J). It was noticed that (j3J) has 
integer solutions only if all three numbers a/ mf + nf, i — 1, 2, 3, have the same irrationality, 
i.e. can be presented as 

.2 i ^2 



mi + nf = ji^/q (4) 

with some integer ^ called weight and the same square-free q called index. In this way, 
set of all wave vectors can be divided into non-intersecting classes Cl q due to class index and 
solutions are to be looked for in each class separately, it is important to realize that this is only 
necessary condition for a solution to exist and some classes can be empty. 

As a first step, we compute the set of all possible indexes q. Due to Lagrange theorem 
on presentation of an integer as a sum of two squares we conclude that q should not be not 
divisible by any prime of the form p = 4u + 3 which reduces the full search substantially. 
Special algorithms for representing square-free numbers as sums of two squares are known, and 
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one of them [18] was used in numerical implementation of our algorithm to compute the set 
of all possible numbers qjf such that all g. t satisfy weight equation (JHJ). Special number- 
theoretical considerations allowed us to disregard a lot of classes from computations (about 
74% of all classes in the domain m, n < 1000) as being empty. 
Next we have to find integer solutions of the weight equation 

111 

- + - = -• (5) 

7i 72 73 

At this step number of variables is reduced from 6 to 3; their individual degrees from 2 to 1, 
and we got rid from irrationality in (j^J). Solutions are looked for only for indexes found at the 
previous step. At this step all solutions of the first equation of © are already found. Finally, 
we check linear conditions on n*. 

In MATHEMATICA implementation, standard functions for list operations and some number- 
theoretical function, like SquareFreeQ and SumOf SquaresRepresentations, from the standard 
package " Number Theory Functions" are used (for details see [19J). 



3 Naive graph presentation 

The graphical way to present 2D-wave resonances suggested in [7] for 3-wave interactions is to 
regard each 2D-vector as a node of integer lattice in spectral space and connect those nodes 
which construct one solution (triad, quartet, etc.) We demonstrate the result in FigH] at the 
upper panel. Obviously, geometrical structure is too nebulous to be useful even in relatively 
small spectral domains. On the other hand, topological structure shown in Figfj] (lower panel) 
is quite clear and gives us immediate information about dynamical equations covering behavior 
of each wave cluster. 

Indeed, energy transport is covered by standard dynamical system, written for simplicity 
for real-valued amplitudes, 

A 1 = a x A 2 A 3} A 2 = a 2 A x A 3 , A 3 = a 3 A x A 2 (6) 

in case of a "triangle" group called further a primary element: {A\, A 2 , A 3 ); by 

A x = a x A 2 A 3 , A 2 = a 2 A x A 3 , A 3 = -{a 3 A x A 2 + a 4 A 5 A 6 ), A 5 = a 5 A 3 A 6 , i 6 = a 6 A 3 A 5 (7) 

in case of "butterfly" group (two connected triangle groups): (A x , A 2 , A 3 )(A 3 , A 5 , A e ), and so 
on. All isomorphic graphs presented in FigH] are covered by similar dynamical systems, only 
magnitudes of interaction coefficients a, vary. However, in general case thus defined graph 
structure does not present dynamical system unambiguously. Consider Figj2] below where two 
objects are isomorphic as graphs. However, the first object represents 4 connected primary 
elements with dynamical system 

{A U A 2 ,A 3 ), (A h A 2 ,A 5 ), {A U A 3 ,A A ), {A 2 ,A 3 ,A 6 ) (8) 

while the second - 3 connected primary elements with dynamical system 

{At,A 2 ,A 5 ), (A^A^A^), (A 2 ,A 3 ,A 6 ). (9) 

To discern between these two cases we set a placeholder inside triangle not representing 
a resonance, we call it further empty 3-cycle. This means that to determine isomorphism of 
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Figure 1: Upper panel: Example of geometrical structure, spectral domain D = 50. Lower 
panel: Example of topological structure, spectral domain \k{\ < 50. The number in brackets 
shows how many times corresponding cluster appears in the chosen spectral domain. 




Figure 2: Example of isomorphic graphs and unisomorphic dynamical systems 



dynamical systems we have to regard graph G together with some parameter(s) 7 to identify 
corresponding dynamical system uniquely. We call a pair (G, 7) an i-pair if it provides isomor- 
phism of dynamical systems. The set of possible parameters 7 (not exhaustive, of course) is: 
number of vertices, their multiplicities, number of edges, their multiplicities, number of primary 
elements (non-empty 3-cycles) N, the list of non-empty 3-cycles L c , etc. Some preliminary 
study of the parameter set show that L c and N is a good first choice, providing a balance 
between informativeness and complexity of numerical implementation. 

Consider a structure (G t , L c ) consisting of: 

• a graph G t each edge of which belongs to at least one 3-cycle; 

• nonempty list L c of some 3-cycles of length 3 of G t such that each edge of G t belongs to 
some cycle(s) of L c . 

Notice that L c does not contain "wrong" triangles. Notation G t has been chosen in order to 



point out that our graphs are, to say, "triangle" graphs. 



Def.l The number of elements in L c is called the order Gt and denoted as N(Gt). 3-cycles 
of G t not belonging to L c are called empty cycles. Number of occurrences of each vertex 
v G Gt in L c is called vertex multiplicity and denoted as n{v). Number of different vertices 
in L c is denoted as M{G t ). 

Obviously, G^-graphs consisting of 3-cycles only, have a very special structure. Our idea is 
to construct a set of all possible graphs of this type of order < iV for some given N inductively, 
beginning with a single triangle. As a next step we can chose all unisomorphic graphs from 
this set and compare corresponding lists L c to find all different dynamical systems. 

3.1 Triangle gluing 

The possibilities of gluing a new triangle to a Gt -graph are not numerous and can be classified 
as follows. Let the new iV-th triangle be T — {v%, i>2, v%\. 

• Vertex gluing. In this case, 1, 2 or 3 vertices of the new triangle are identified with (glued 
to) vertices of some distinct triangles of the graph, constructed at previous inductive step, 

• Edge gluing. In this case, 1, 2 or 3 edges and corresponding vertices of the new triangle 
glued to edges and vertices of some distinct adjacent triangles of the graph. Notice that 
gluing of three edges is simply filling an empty triangle of G?j; ~ . 

• Mixed gluing. In this case, one vertex vn 1 of the new triangle is glued to a vertices of 
some triangle of the graph and the edge vn 2 vn 3 is glued to an edge of another triangle. 

The cases described above are illustrated by figures FigJS]- [H 




a) Gluing by one vertex 




b) Gluing by two vertices 




c) Gluing by three vertices 



Figure 3: Vertex gluing of a new triangle to G t 
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a) Gluing by one edge 

b) Gluing by two edges 

+ $ * $ 

c) Gluing by three edges 
Figure 4: Edge gluing of a new triangle to G[ N ^ . 

Figure 5: Mixed gluing of a new triangle to G\" l> . 

• By vertex gluing, the G\ N ^ structure is enhanced by: 

— two vertices and three edges (one vertex glued) 

— one vertex and three edges (two vertices glued) 

— three edges (three vertices glued) 

• By edge gluing, the G^ JV ~ 1) structure is enhanced by: 

— one vertex and two edges (one edge glued) 

— one edge (two edges glued) 

— the graph structure stays unchanged (three edges glued) 

• By mixed gluing, the G[ N ~^ structure is enhanced by two edges. 

In each case, the list of graph G^' 1 ^ is extended by the vn 1 vn 2 v n 3 cycle of the new 

triangle (or whatever these vertices will be called after the gluing). 

3.2 Some estimations 

Enhancing a given G t with V vertices and E edges by a triangle, we encounter the following 
possibilities: 
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• a sole triangle not connected to the existing graph is added (1 possibility); 

• vertex gluing (V possibilities); 

• edge gluing (E possibilities); 

• mixed gluing - approximately E/V(E/V — 1) possibilities; 

• filling an empty triangle (very rare). 

Therefore at each inductive step the mean number of vertices is V < 1.5N and the number of 
arcs can be roughly estimated as E < 2N. Therefore, the number of emerging unisomorphic 
graphs can be estimated from above as some ~ 4N and the overall number of graphs at step 
N is 0(N 2 ). 

4 Hypergraph presentation 

To diminish computational time and complexity we construct a hypergraph presentation of i- 
pairs introduced in the previous section. A hypergraph is a structure that consists of a set of 
vertices and a multiset of edges, called hyperedges. A hyperedge is a set of vertices, all vertices 
in such a set are connected. The collection of hyperedges is a multiset because it is possible that 
some hyperedges appear more times. A traditional graph is a special case of a hypergraph, in 
which all edges are two-element sets and do not appear more than once. For the representation 
of 3-wave resonances we consider the triangles as "the nodes" of corresponding hypergraph. 

Def. 2 A hypergraph with 3-cycles of a triangle graph Gt as its vertices and nodes (m, n) 
of Gt as its edges is called a triangle hypergraph and is denoted as HGt- Set of its vertices 
and edges is denoted as Vug an d Ehg correspondingly, i.e. HG t = (Vug, Ehg)- 

Notice that since a node (m, n) of G t can belong to several 3-cycles, corresponding HG t 
has in fact hyperedges instead of edges of a simple graph. A hypergraph HG t generated by 
Gt has two properties: 

• Each vertex is part of exactly three hyperedges. 

• Each pair of vertices is part of at most two hyperedges. 

The first property follows from the fact, that each vertex of HGt represents a 3-cycle which 
consists of three different nodes of G t - If the second property is violated then the two associ- 
ated 3-cycles of G t have three nodes in common, hence they are identical. 

As an illustrative example, let us write out explicitly a hypergraph presentation of the 
dynamical systems (JE]) and ([9]) presented in Figj2] at the left and right panel correspondingly: 

(V HG = {1,2, 3, 4}, E HG = {{2}, {3}, {4}, {1, 2, 3}, {1, 2, 4}, {1, 3, 4}}) (10) 

and 

(V HG = {1,2, 3}, E HG = {{1}, {2}, {3}, {1, 2}, {1, 3}, {2, 3}}) . (11) 
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4.1 Incidence matrix 



For computation purposes it is convenient to represent a hypergraph HGt by its incidence 
matrix which is constructed in the following way. 

Def. 3 A rectangular matrix J = (fij) with M(G t ) columns and N(G t ) rows is called 
incidence matrix of G t if 



fi,j 



1, j — th non-empty 3-cycle contains i — th node, 
otherwise. 



(12) 



Each column of the matrix $ represents a triangle in the solution set of (E]) while each row 
represents a node (see Def. 1). Since we are not interested in nodes themselves but in their 
relation to each other we can relabel the nodes of the triangle with ascending integers in an 
arbitrary way and use the labels of the nodes for indexing elements in a matrix. Now we can 
construct the hyperedges of HG t : if the j-th entry of a row is equal to 1 then we add j 
to this hyperedge. The vertices of HGt are elements of L c . The ordering of the hyperedges 
is not important, because it is a multiset. However, it is better to have a "normal form", so 
we sort the hyperedges by using some ordering. Since we are interested in an implementation 
in MATHEMATICA we choose the ordering used by the command Sort. This is an ordering, 
which orders lists ascending by their length, and lists of same length lexicographical by their 
elements. For the dynamical systems there is no ordering with practical advantages for the 
implementation, so we let them unsorted. The incidence matrices of dynamical systems (IE]) 
and ([9]) have form 



/ 1 
1 
1 





i 
i 




1 





o\ 
1 
1 




1 



(13) 



and 



/ 1 
1 


1 



o\ 

i 

i 




1 J 



(14) 



correspondingly. Analogously, incidence matrices of their hypergraphs, here we use the ordering 
of the hyperedges described above, 



/ 1 \ 



1 

1 1 




1 




1 

110 1 
\1 1 1/ 



(15) 
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and 



/ 1 


o 


\ 


o 


1 


o 


o 


o 


1 


1 


1 


o 


1 





1 


Vo 


1 


1 / 



(16) 



are also different. The incidence matrices of the dynamical system and the corresponding hy- 
pergraph are not identical. The reason is the use of different orderings for the vertices. Matrix 
( fTo*]) is just a permuted version of matrix (fTB"]) . Since we use a special ordering for the hyper- 
edges, which are described by the rows of the incidence matrix, we obtain permuted rows. For 
identifying isomorphic dynamical systems it is not necessary to preserve an ordering, because 
dynamical systems with permuted elements are still isomorphic. Hence, neither row permuta- 
tions nor column permutations destroy the isomorphism of dynamical systems. In this example 
only row permutations occur, permutations of columns are just another ordering of the elements 
of L c . This construction can be redone and the dynamical system can be reconstructed out 
its hypergraph: by considering the columns of this matrix we know which nodes belongs to a 
certain 3-cycle. 

Obviously, if two hypergraphs ([TO]) and fTTTl) are not isomorphic, also their incidence ma- 
trices ([TBI and ff"Hj) are different. But in general for the final decision it is necessary to have 
an algorithm to establish isomorphism of hypergraphs. Since there are not so many general 
algorithms for hypergraphs one has to find a representation where it would be possible to use 
standard algorithms for graph isomorphism. This leads us to auxiliary multigraph construction 
presented in the next section. 

4.2 Multigraph construction 

A multigraph MGt is constructed in the following way. Its vertices coincide with the vertices 
of HGt and each hyperedge is replaced by all two-element subsets. To maintain the whole 
information we have to label the created edges so that edges which belong to the same hyperedge 
of HG t are labeled identically. These labels allow to reconstruct HG t and J which is a 
necessary step while generating dynamical systems. The hyperedges which contain only one 
vertex can be omitted because they contain no further information about the cluster structure. 
Of course, some edges may occur in MG t twice - this is the case if two 3-cycles of G t share 
two nodes. Figure [6] shows two multigraphs corresponding to the dynamical systems shown in 
FigfSJ For easier distinction we use triangle symbols for the vertices of the multigraphs, because 
a vertex represents a 3-cycle of G t . 

A multigraph MGt has following properties: 

• At most two edges connect a pair of vertices. 

It follows from the fact that a pair of 3-cycles can share at most two nodes. If they would 
share also their third node, they would be identical. 

• At most three different labeled edges can occur at a vertex. 

A 3-cycle has three nodes therefore it can only share three different nodes with other 
3-cycles. 

• The number of vertices is equal to the number of non-empty 3-cycles in G t . 
By definition. 



10 




Figure 6: Multigraph representations for dynamical systems (JEJ) and correspondingly 

• The total number of edges with identical labels is ( - p ~ 1 ' )p , where p is the number of 
elements in the corresponding hyperedge. 

Edges with identical labels belong to the same hyperedge, and the number of two- element 
subsets is ( p ) . 

4.3 Hypergraph versus naive graph 

Summarizing briefly the procedure described above, following has been done: 

• all integer solutions of (EJ) are found; 

• topological presentation of the solution set as an i-pair (G t , L c ) is constructed which 
presents corresponding dynamical system uniquely up to isomorphism; 

• i-pair (G t , L c ) is transformed uniquely into a hypergraph HG t ; 

• for computational purposes, some auxiliary builds are introduced - incidence matrix 
T{G t ) and multigraph MG t ; both maintain the isomorphism of dynamical systems. 

The advantages of hypergraph representation compared with a more simple i-pair represen- 
tation given in Sec. 2 are following: 1) no additional parameter to distinguish non isomorphic 
dynamical systems are needed; 2) a standard graph isomorphism algorithm can be used to es- 
tablish the isomorphism of multigraphs; 3) the size of constructed multigraphs is approximately 
one half of that for Gt- Some results of MATHEMATICA implementation of this procedure 
are given in the next Section. 

5 MATHEMATICA implementation 

Details of our MATHEMATICA implementation can be found in [19] (solutions of (CD) and 
geometrical structure) and in [20] (topological structure and dynamical systems). General 
computation schema is following. We implemented algorithm sketched in SecJ21 computed 
all solutions of CQ) and used MATHEMATICA package "DiscreteMath'Combinatorica"' to 
plot triangle graph Gt, and to construct incidence matrix $ and multigraph MGt- To 
establish multigraph isomorphism we modified a standard algorithm provided by the "Dis- 
creteMath'Combinatorica"' package, because it can only be used for for simple graphs and 
multigraphs with unlabeled edges. Some necessary conditions of multigraphs isomorphism are 
checked as a preliminary step, in order to make computations faster. As an output, list of all 
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resulting clusters is given, for each of them corresponding incidence matrixes, hypergraphs, and 
dynamical systems are written out, and graphs G t and MG t are plotted. We also compute 
how many isomorphic clusters of each form appear in the chosen computation domain. Results 
for computation domain D = 50 are given below. 

1. 18 systems: ( V HG = {1}, E HG = {{1}, {1}, {1}} ) 

1 \ A x = aiA 2 A 3 
1 A 2 = a 2 A 1 A 3 

1 / A 3 = a 3 A x A 2 



A 




2. 4 systems: 



/ 1 
1 



V 1 





1 
1 

1 / 



V HG = {1,2},E HG ~- 

A x = a x A 2 A 5 
A 2 = a 2 A x A^ 
A 3 = a 4 A A A 5 
A A = a 5 A 3 A 5 
A 5 = \ {a 3 A x A 2 + a 6 A 3 A 4 



{!},{!}, {2}, {2}, {1,2} 





3. 1 system: 



/ 1 
1 




1 






1 
1 



1 / 



V HG = {1,2,3}, E HG ~- 

Ai = a l A 2 A 6 
A 2 = a 2 AiA 6 
A 3 = a 4 A 6 A 7 
A 4 = a 7 A 5 A 7 
A 5 = a 8 A 4 A 7 

A 6 = | {a 3 A x A 2 + a 5 A 3 A 7 ) 
A 7 = \ (a 9 A 4 A 5 + a 6 A 3 A 6 ) 



{!},{!}, {2}, {3}, {3}, {1,2}, {2, 3} 
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2 systems: 


v u 

V 1 


/ 1 o 





o \ 


1 








1 








1 











1 











1 








1 





1 


1 


V 1 1 


1 


/ 



HG 



= {1,2,3,4}, E HG ~- 

A\ = a\A 2 A 9 

A 2 = a 2 A x A 9 

A 3 = a 4 A 4 A 9 

At = a 5 A 3 A 9 

A 5 = a 7 A$A 9 

A 6 = a 10 A 7 A 8 

A 7 = a u A 6 A 8 

A s = | (a 12 A 6 A 7 + a 8 A 5 A 9 ) 

^9=1 {a 3 AiA 2 + a & A 3 A± + a 9 A 5 A 8 ) 



{1}, {1}, {2}, {2}, {3}, {4}, {4}, {3, 4}, {1, 2, 3} 




1 syst 


em: ^ 


V HG 


/ 1 








o\ 




1 
















1 
















1 
















1 













1 




1 


1 













1 


1 







\o 





1 


1 J 





{1,2,3, 4}, E HG ~- 

A\ = a\A 2 A 7 

A 2 = a 2 A x A 7 

A 3 = (x±A 7 A 8 

A 4 = a 7 A 8 A 9 

A 5 = a w A e A 9 

A 6 = a u A 5 A 9 

A 7 = \ {a 3 A x A 2 + a 5 A 3 A 8 

A 8 = | {a 6 A 3 A 7 + a 8 A 4 A 9 

A 9 = ~ (a 12 A 5 A 6 + a 9 A A A 



{1}, {1}, {2}, {3}, {4}, {4}, {1, 2}, {2, 3}, {3, 4} 
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6. lsystem:( V HG = {1, 2, 3, 4, 5}, E HG = {{1}, {1}, {2, 4}, {2, 5}, 
{3,4}, {3,5}, {4,5}, {1,2,3}} ) 

Ai = aiA 2 A 8 
a 2 A 1 A s 

| (ai v4 5 v4 7 + a 4 v4 4 v4 8 ) 

■ (ai 3 A 6 A 7 + a 5 A 3 A 8 ) 

■ (a u A 3 A 7 + a 7 A 6 A 8 ) 
_ (a u A A A 7 + a 8 A 5 A 8 ) 

A 7 = \ {a 12 A 3 A 5 + a 15 A 4 A 6 ) 
i 8 = | (a 3 AiA 2 + a 6 A 3 A 4 + a 9 A 5 A 6 ) 



( 1 
1 

1 
1 




V i 



\ 


1 









1 
1 



1 1 





1 



1 
1 

o ) 



A 2 
A, 
A, 
A, 
A, 



2 
1 
2 

_ 1 




7. 1 system: I Vj 



{2}, {3}, {4}, {7}, {7}, 



/o 


1 














o\ 








1 
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1 

















1 





1 

















1 











1 
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1 
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1 
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At = a 4 A 6 A 8 
A 2 = a 7 A 7 A 12 
A 3 = a 10 A 9 A 12 
A 4 = a 19 A 5 A n 
A 5 = a 29 A A Au 
A & = \ (a 5 AiA 8 + axA 7 A 12 ) 
A 7 = \ (a 8 A 2 A 12 + a 2 A 6 A 12 ) 
A 8 = \ {a§A x A % + a ie A 10 A n ) 
A 9 = \ (a u A 3 A 12 + a 13 A w A 12 ) 
A w = \ (a 17 A 8 Au + a 1A A 9 Ai 2 ) 
An = \ (a 2 iA 4 A 5 + ai 8 A 8 A w ) 
An = \ (a 9 A 2 A 7 + a 3 A e A 7 + a u A 3 A 9 + ai 5 A 9 A w ) 
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These results show that in spectral domain D = 50 which contains ~ 2 • 10 3 Fourier 
harmonics we have only 7 non-isomorphic dynamical systems (clusters of waves) for further 
analytical and numerical study. Some of them, for instance ([6]), are known to be solved explicitly 
in Jacobean elliptic functions (example of explicit expressions for the case of spherical planetary 
waves can be found in [21] )• Knowledge of the explicit form of a dynamical system allows 
sometimes to obtain a few conservation laws as in the case ([7]) and simplify substantially further 
numerical investigations of these systems. It is important to understand that though qualitative 
properties of all isomorphic clusters are the same, their quantitative properties depend on the 
magnitudes of coupling coefficients ati, of course. Computation of these coefficients is usually 
done by standard multi-scale method which is tedious but completely algorithmic procedure 
and can also be programmed in MATHEMATICA (see [TH] for its implementation). 




6 Two mechanisms to destroy clusters 

There are two mechanisms which can destroy clusters constructed above: 1) increasing of 
spectral domain D, and 2) taking into account quasi-resonances, i.e. integer solutions of 

01±0 2 ±...±0 s = ft >0, fci±jfe 2 ±....±ife« = (17) 

with some non-zero resonance width Q. Below we regard briefly both of them. 



6.1 Increasing of spectral domain 

Obviously, the structure of clusters becomes simpler with diminishing of the domain D - some 
solutions (triads) disappear. On the other hand, increasing of D might lead to a substantial 
changes of the structure. Thus it is important to understand how solution structure depends 
on the chosen computation domain. With this aim let us re-write first equation of ([3]) in the 
form 

111 

s + srs (18) 

and notice that k 3 < k\ and k 3 < k 2 . Introducing notations k_,k ,k + for the minimal, 
intermediate and maximal of the numbers k±,k 2 , k% we conclude that k^ = k^. Without loss 
of generality one can assume that k 2 < k\ and let k 2 = &o and k\ = k + (it is a formal 
notation taken for convenience of calculations, because k 2 and k\ can also coincide). Now 
let us fix /c_ and re-write (fTSl) as 

k k k_ 
+ k — k_ 

Last expression achieves maximum if ko = k^ + 1 which yields k + < k^ik^ + 1) and similar 
considerations show that also ko < k^(k^ + 1). This means that wave interactions are local 
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in the following sense: lengths of wave vectors constructing a solution of (EJ) can not be too 
far apart. In particularly, if we are interested in the solutions structure in the domain, say, 
ki < D = 50, it is enough to investigate a larger domain D = 50 2 + 50 = 300, in order 
to establish which clusters stay unchanged and to find those which are enhanced via solutions 
with wave vectors lying outside of the initial domain D = 50. 

6.2 Quasi-resonances 

It was shown in [23] that for discrete quasi-resonances to be able to start some low boundary 
for resonance width Q can be written out explicitly. It interesting that for many dispersion 
functions, there exist global low boundary for most clusters which does not depend on the 
spectral domain under consideration and also does not depend on the number of interacting 
waves s. For instance, in case of = (m 2 + n 2 ) 1//4 (gravity water waves) the use of the 
generalized Thue-Siegel-Roth theorem [21] yields Q > 1. Obviously, for arbitrary dispersion 
function a local low boundary exists which is defined by the spectral domain T = {(m, n) : < 
\m\, \n\ < D < oo} chosen for numerical simulations. Indeed, let us define = min p Q p , 
where 

n p = \ (k p 1 )±0(k p 2 )±... ±0(#)|, $=(TnJ,nJ)eT, Vj = l,2,..., S , 

and 

0(k p 1 )±0(k p 2 )±...±0(k p s )^O Vp, 

and index p runs over all wave vectors in T , i.e. p < AD 2 . So defined Q p obviously is 
a non-zero number as a minimum of finite number of non-zero numbers and Qd is minimal 
resonance width which allows discrete quasi-resonances to start, for chosen D. 

Physically important resonance width Q p h ys is defined by the accuracy of computations and 
precision of measurements in numerical and laboratory experiments correspondingly. Quasi- 
resonances with Q D > Q p hys will not destroy the clusters. 

7 Discussion 

In order to apply theory of discrete resonances to a real physical problem, the profound study of 
constructed dynamical systems is needed. It is well-known that dynamical system demon- 
strates periodic energy exchange between the modes of a triad. On the other hand, dynamical 
systems consisting of a few connected triads have enough degrees of freedom M to behave 
chaotically. The question of major importance therefore is to discern between two classes of 
situations: 1) resonance clusters with periodic energy exchange within each cluster, and 2) 
those which can be described statistically, similar to kinetic equation approach. From this 
point of view, all our theoretical results and symbolical programming can be regarded as an 
introductory step for further numerical simulations. 

We are quite aware of the fact that there exists multitudinous number of important questions 
to be answered in order to understand a very complicated mutual relationship between discrete 
and statistical regimes of wave system dynamics. For instance, is corresponding statistical 
dynamics close to Gaussian? Is a probability of attractor appearance in the subspace generated 
by integrals of motion uniformly distributed? What is the minimal value of M allowing to 
"forget" topological details of the discrete, low-dimensional dynamical system and describe 
corresponding dynamical system statistically? How does energy exchange between isolated and 
continuous subsystems look like? What is the role of nonlinearity in triggering energy flux 
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toward small scales? Is it possible to develop some analytical tools for description of low- 
dimensional systems (for example, generalized kinetic equation that accounts the finite width 
of frequency resonances)? etc. 

A first feeling of possible answers to some of these questions can be obtained by computer 
simulations with a few well-chosen dynamical systems with degrees of freedom from TV = 4 
to 10 4- 20 which is in our agenda. Notice that in the case of 3-wave resonances one has 
to construct dynamical systems as it was done above and choose which are not enhanced by 
increasing of computation domain. Choice of initial conditions for numerical simulations would 
be another important subject to study for, as it was mentioned in [21], even for one isolated 
resonant triad it is always possible to chose initial energy distribution among the modes in such 
a way that the period of their energy exchange will tend to infinity. 

In general our graph-theoretical approach can be used, with appropriate refinements, also for 
s-wave resonances, with s > 4. In this last case, also existence of different types of resonances, 
spectrum anisotropy, etc. [23] have to be taken into account in order to choose representative 
dynamical systems for numerical simulations. 

The same approach (algorithms from [TS], graph construction, etc.) can also be used directly 
for any mesoscopic system with resonances of a more general form 

Pi0i ±p 2 02 ± ••• ±p s 0s = 0, pih ±p 2 k 2 ± .... ±p s k s = (19) 

with integer pi 

Acknowledgements. Authors acknowledge the support of the Austrian Science Founda- 
tion (FWF) under projects SFB F013/F1304 and SFB F013/F1301. 

References 



[2] 
[3] 
[4] 
[5] 
[6] 

[7] 
[8 
[9 
[10 

[11 
[12 



P. Denissenko, S. Lukaschuk, S. Nazarenko. E-print: arXiv.org:nlin. CD/0610015 (2006) 

R.Toral, C.J.Tessone. Comm. Comp. Phys. 2 (2), pp. 177-195 (2007) 

A. Stefanovska, M.B. Lotric, S. Strle, H. Haken. Physiol. Meas. 22, pp.535-550 (2001) 



H. Spohn. E-print arXiv:math-ph/ 0605069 (2006) 



E.A. Kartashova. Phys. Rev. Let. 72, pp. 2013-2016 (1994) 

V.E. Zakharov, V.S. L'vov, G. Falkovich. Kolmogorov Spectra of Turbulence, Series 
in Nonlinear Dynamics (Springer, 1992) 

E.A. Kartashova. AMS Transl. 2, pp. 95-129 (1998) 

C. Connaughton, S. Nazarenko, A. Pushkarev. Phys. Rev. E 63, 046306 , 4 (2001) 

Yu.V. Lvov, S. Nazarenko, B. Pokorni. Physica D 218 pp. 24-35 (2006) 

M. Tanaka. J. Fluid Mech. 444, pp. 199-221 (2001) 

M. Tanaka. Fluid Dynamics Research 82(1), pp. 41-60 (2001) 

M. Tanaka, N. Yokoyama. Fluid Dynamics Research 34(3), 199-216 (2004) 



17 



[13] M. Tanaka. To appear in J. Phys. Oceanogr. (2007) 

[14] V.E. Zakharov, A.O. Korotkevich, A.N. Pushkarev, A.I. Dyachenko. JE TP Letters, 82 (8), 
pp. 487-491 (2005) 

[15] E. Kartashova. Letters JETP 83 (7), pp. 341-345 (2006) 

[16] E. Kartashova. J. Low Temp. Phys. 145 (1-4), pp. 286-295 (2006); E. Kartashova, A. 
Kartashov: Int. J. Mod. Phys. C 17(11), p p. 1579-15 96 (2006); Comm. Comp. Phys. 2 
(4), pp. 783-794 (2007); E-print: arXiv.org jmath-ph/0701030[ To appear in Physics A: 
Stat. Mech. Appl. (2007) 

[17] E.A. Kartashova, G.M. Reznik. J. Oceanology 31, pp. 385-389 (1992) 

[18] J.M. Basilla. Proc. Japan. Acad. 80, Ser. A, 40 (2004) 

[19] Ch. Feurer, E. Kartashova, G. Mayrhofer, C. Raab, W. Schreiner. RISC Report Series. 
March, 2007 

[20] E. Kartashova, G. Mayrhofer. RISC Report Series. March, 2007 



[21] E. Kartashova, V.S. L'vov. E-print arXiv.org.nlin/0606058, Submitted (a shortened ver- 
sion) to Phys. Rev. Let. (2007) 

[22] S. Skiena. Implementing Discrete Mathematics. Addison- Wesley, Redwood City, CA 
(1990) 

[23] E. Kartashova. E-print arXiv.org:math-ph/070107. Submitted to Phys. Rev. Let. (2007) 

[24] W.M. Schmidt. Diophantine approximations. Math. Lecture Notes 785, Springer, 
Berlin, 1980 



18 



